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Abstract: The EIT waves [named after the Extreme-ultraviolet Imaging Telescope 
(EIT) onboard the Solar and Heliospheric Observatory (SOHO)] are in the literature 
usually described as fast magneto-acoustic (FMA) modes. However, observations show 
that a large percentage of these events propagate with very slow speeds that may be as 
low as 20 km/s. This is far below the FMA wave speed which cannot be below the sound 
speed, the latter being typically larger than 10^ km/s in the corona. In the present study 
it is shown that, to account for such low propagation speed, a different wave model should 
be used, based on the theory of gravity waves, both internal (IG) and surface (SG) ones. 
The gravity modes are physically completely different from the FMA mode, as they are 
essentially dispersive and in addition the IG wave is a transverse mode. Both the IG and 
the SG mode separately can provide proper propagation velocities in the whole low speed 
range. 



PACS Numbers: 96.60.P-; 52.30.Ex; 52.35.Fp; 



I. INTRODUCTION 



One of the most striking details related to so called EIT waves is the wide range of the 
propagation velocity reported in the literature. Velocities of 50-60 km/s were reported 
in Ref. 1, together with those of 350 km/s. In the same reference, values of 500 km/s 
are mentioned also as observed by Yohkoh/SXT. In Ref. 2 the values in the range 25-450 
km/s are reported observed by SOHO EIT. Much more on that issue is available in Ref. 3, 
and also in Ref. 4 where velocities as low as 20 km/s are reported. 

The EIT wave is usually described in terms of the fast magneto-acoustic mode, which 
implies the phase propagation speed given by 



Clearly, the Alfven velocity Ca, determined by two plasma parameters gives some space 
to adjust to such a large span of values, and even to provide the low values mentioned 
above. For example, taking uq = 10^^ m^^/^, Bq = 10~^ T, and ng = 3 ■ 10^^ m"^'''^, 
Bq = 2 - 10^^ T, for Ca yields 22 km/s, and 25 km/s, respectively. It remains questionable 
though how realistic such parameters are. Even the 'classic' value of the phase speed of 
the EIT waves, that is around 250 km/s, implies tiq = 2 ■ 10^^ m~^^^, Bq = 2 ■ 10~^ T, 
hence, a rather low value for the magnetic field. The impression is that, to account for 
such speeds, certain values of the plasma-/? are just adopted with the only purpose of 
fitting it into the assumed and pre-selected fast magneto-acoustic wave model. However, 
according to Ref. 3, out of 160 studied events of the EIT waves, more than 60 percent 
were with the speed below Ca- 

On the other hand, the FMA wave speed includes the sound speed too. In calculating 
Cs one has much less freedom for adjustment, because it is determined by the temperature 
only. For example, taking the ion and electron temperatures equal, and changing it in the 
range 0.6 — 2 million K, one obtains in the range of 100 — 182 km/s. Thus, remains 
above 10^ km/s, and no adjustment of the plasma-/? may help. According to Eq. ([T]) there 
will always be Vp > Cs, therefore the observed low values of Vp mentioned above (just a 
few tens of km/s) can not possibly be explained within the fast magneto-acoustic wave 
model. 

This inspired some researchers to make a completely different models, claiming that 
the EIT perturbations are not waves at all. Such non-wave models are merely focused 
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on how the EIT perturbations are produced in the first place, e.g., due to some specific 
change in the magnetic field topology. However, the origin of perturbations is not an 
issue here. Prom basic plasma theory we know that an initial perturbation, no matter 
how produced, will propagate as a wave, and regardless of its shape, it can be decomposed 
into Fourier harmonics. Hence, one has to deal with waves in any case. 

What remains as one option is to search for some different wave model that could give 
proper propagation speed values in the low velocity range mentioned above. This should 
by no means reduce the validity of so widely accepted fast magneto-acoustic (FMA) 
wave model, which indeed seems to be able to capture many essential features of the 
fast EIT perturbations. Rather, it should represent a possible supplement in the slow 
speed domains where the FMA model appears to be not applicable. Regarding the waves 
propagating orthogonal to the magnetic field, that seems to be the case with the EIT 
perturbations, the choice is very limited, especially in view of the time and space scales 
involved in the perturbations. Yet, such plasma modes exist, and those are the internal 
gravity (IG) wave and the intermediate surface gravity (SG) wave, which have features 
completely different as compared to the FMA mode. Those features of interest for the 
EIT waves will briefly be described in the following sections. 

II. INTERNAL GRAVITY WAVES 

The physics of the IG wave is well known from standard books dealing with fluids 
in external gravity fleld, as excellent examples to mention just the books of Kundu^ and 
Sturrock^, where the mode is described in the frame of ordinary single fluid theory. In 
the present two-component plasma case applicable to the coronal plasma, the starting 
equations for ions, in usual notation, are the continuity, momentum and energy equation, 
respectively, 

^ + V • {nSi) = 0, (2) 
miUi I ^ • V j iTi = -eniVcf) - Vpi + miUig, (3) 



(4) 



Here, the gravity is assumed in the 2;-direction, g = —gCz, the ions are singly- charged, 
and is the electrostatic potential. 
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For low-frequency perturbations (with a phase speed much below the electron thermal 
velocity), electrons will be treated as isothermal, and their motion will be described by the 
quasi-neutrality condition ng = rti (implying spatial scales far above the Debye length), 
and by the momentum equation in the inertia-less limit 

= engVcf) — KTeVne + meUeg- (5) 

The magnetic field is neglected and a proper justification for this will be given in the 

forthcoming text. 

A. Equilibrium plasma 

Assuming a static and isothermal equilibrium without electric fields, and for a constant 
gravity, for the density from Eqs. ([3], one obtains the barometric formula 

Hj exp{—z/H), (6) 

where H = KTj/{mjg) {j = e, i) is the characteristic scale length for the stratified medium. 
Note that the same expression can be obtained for an ordinary neutral gas. According 
to the formula (jH]), there should be a difference (with height) in the density for species 
with different mass. In the case of electrons and ions this would then violate the quasi- 
neutrality. To preserve the quasi-neutrality, the plasma will produce a macroscopic elec- 
tric field Eq = EqCz [known also as the Pannekoek-Rosseland (PR) field^'* or ambipolar 
gravity-induced field] which opposes the gravity in order to maintain its macroscopic quasi- 
neutrality. The electric field is oriented in such a way to lift the ions against the gravity, 
and in the same time to pull the electrons down. Hence, the appropriate momentum 
equations for the two species will read 

eno-Eo + hTq^^ + meUog = 0, Eq = -V0o, (7) 
dz 

enoEo - kTq^ - miUog = 0. (8) 
dz 

Here, we used the quasi-neutrality n^o = Ueo = uq, and for simplicity the same temperature 

for both species is assumed. Subtracting the two equations yields the PR field 

^ (nii — me)q rriiq , , 

Eq = ^ ~ 9 

° 2e 2e ^ ' 

Setting this into Eqs. ( ITfHj) . contrary to Eq. ([6]) one now finds the same distribution for 

both electrons and ions 

no = iVoexp(-^), = (10) 



It is seen that the plasma characteristic length Hp is twice larger than H for a neutral gas 
(obtained above also for ions without the effects of the PR field), and in the same time 
for electrons it is shorter than if by a factor 2me/mj. These conclusions are general and 
valid for any external gravity field and plasmas with two species. Yet, these expressions 
are approximative because g and Tq are assumed locally constant with height z. An 
application of this effect to ion acoustic oscillations can be found in Refs. 9, 10. It can 
easily be shown that for two general species a and b, with the charge —Qa and qt,, and 
with different temperatures 7^ Tj,, the PR field will read 

rribTa-maTb 

Eo = g — ■ (11) 

QaTb + QbTa 

Although frequently overlooked in the literature, this electric field is reality. The field of 
the same origin has been experimentally detected on the surface of metallic conductors. ^^'^^ 
It appears due to the redistribution of electrons and nuclei in metallic objects placed in 
the gravity field. In the case of an expanding atmosphere (e.g., in the solar wind) the 
magnitude of the electric field can substantially be greater than the Pannekoek-Rosseland 
value presented here. Details on that issue are available in Refs. 13, 14. 

The presence of several charged ion species may drastically change the previously ob- 
tained monotonous exponential distribution (fTOj) . To demonstrate this, one may assume 
the presence of / singly charged ion species of equal temperature and write the corre- 
sponding equilibrium momentum equations —nTdrij/dz — rijimjg + erijE = 0, j = 1, ■ ■ ■ /. 
Performing a summation of the set of / equations, together with the electron equation, 
with the help of quasi- neutrality, ng = Y^irij, for the ambipolar electric field we obtain 
E = m{z)g/2e. Here, 

/ I 
1 1 

is an effective, spatially dependent mean ion mass. Using these expressions in the equation 
for one specific j-th ion species, we obtain d{lnnj)/dz = —[nij — m{z)/2]g/{K,T). It is 
seen that in fact it yields a growth with z as long as rrij < m{z)/2. Hence, lighter ion 
species may have an inverse density distribution in lower layers. It may be concluded that 
in general Hj = Hj{z), and this dependence includes a possible change of the sign too. 
Such distributions may be expected in environments with dominant heavy ions, details 
related to the terrestrial atmosphere are available in Ref. 15. In the corona the situation 
is just the opposite, therefore the effects of extra ion species will not be studied here. The 
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presented picture of the PR field and the change in the scale length (in the case of equal 
temperatures enlarging it by a factor 2 for protons, and in the same time reducing it by 
a factor 1000 for electrons) are based on the model of a quiet, stratified fluid. However, 
turbulence and mixing of different layers, especially in the lower solar atmosphere, may 
effectively reduce its importance. 
B. Small perturbations 

Linearized momentum equations for the two species, with the help of Eqs. (ITl fTUl) . 
become 

mifiQ^^ = -e?ioV0i - Vpii - nigrrii ^ ^. _ ^z, (13) 



at l + T 

= enoV0i — /tTeVrii — riigrrii e^. r 



l + T 



(14) 



Note that the gravity term in the electron equation contains the ion mass. The gravity 
terms in the two equations are equal if r = 1, implying that the electrons and protons 
effectively behave as particles of the same weight (but with different inertia, hence on 
the left-hand side in the electron equation). This is the consequence of the PR field. The 
other two linearized equations are 

dfii 



dt 

dp I -ypo drii 



+ noV ■ Vii + Vii ■ VriQ = 0, 



+ nl {v,i ■ V) (poUo ' 



0. 



(15) 



(16) 



dt uq dt 

To satisfy flux conservation, the perturbations propagating in the (x, 2;)-plane are 
taken of the form^'^'^*^ 



ni = Hi exp 



Pi = pi exp 



z 



2H 

z 



2H 



+ i {kxX + k^z — Lot) 
+ i {kxX + kzZ — out) 



V{x,z)i = V{x,z)i exp 



— + i {k^x + kzZ - ujt) 
Zn 



)i = <pi exp 



2H 



+ i {k^x + kzZ — (jjt) 



(17) 

(18) 
(19) 
(20) 



where pi etc., is used to denote the constant amplitude. Setting these expressions into 
Eqs. ( 1T51ITB]) yields the following equations for the amplitudes ni^pi^Vxi^Vzi- 



k^nTeUi + k^pi - uniiNoVxi = 0, 
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—iuniiNoVzi = 0, 
-iufii + ik^NoV^i + (^k^ - Nqv^i = 

-iwynTini - iupi H — v^i = 0. 

hp 

Here, a = Ti/(Te + %), the number of equations was reduced by eliminating the elec- 
trostatic potential by using the electron momentum equation. Its x-component gave 
iT'il'no = e0i/(KTe), while its 2;-component was identically equal to zero. Due to the shape 
of Eqs. (I17ti20p . and after using Eq. ffTOj) . the common exponential term cancels out. After 
a few steps this yields the dispersion equation 



Ul + Vl {kl + ^ ^2^2^2 ^ Q (21) 



Here, 



^ mj^i ' ^ n{Te + Ti)(Te + -fTi)' 
In i/p the temperatures of the two species are kept different, Uc is the ion acoustic (lA) 
wave cut-off frequency, and Uf, is the buoyancy [Brunt- Vaisala, BV] frequency. 

Eq. fl2ip describes coupled the lA and the internal gravity (IG) waves. The frequencies 
of the two modes are usually well separated, for comparison see a detailed analysis for 
ordinary fluids in Ref. 6. Hence, neglecting the last term in Eq. fl2ip yields the lA mode 



2 2 



+ {kl + kl) vl (23) 



On the other hand, in the low frequency range one may neglect o;^ in Eq. fl2T]) . which 
yields 

2 _ ^bkx'^s _ ^bkx^s i,2 _ 7,2 , 1.2 

In case when u"^ <^ k'^v'l this yields the frequency for the IG mode presented in Ref. 5 for 
incompressible fluids and after using the Boussinesq approximation, 

-\- rv 
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According to conventional criterion of incompressibility, as long as the perturbed velocity 
|fi| is much below the sound velocity in the system, Eq. f l24|) will be valid. In fact, 
obtaining Eq. ( l24l) was the main reason for the preceding text. Some well known and 
interesting properties of the IG mode can easily be seen by analyzing Eq. (!24l) . The phase 
velocity of the mode is 

Via = = ^3 [k^ex + kzGz) ■ (25) 
On the other hand, its group velocity is 

C/G — ^^'q^ ~ ^3 {kzf^x kx^z) ■ (26) 

This yields 

CjG ■ = 0. (27) 

From Eqs. ( I25ll27|) it is seen that the horizontal components of Cig,Vig are in the same 
direction, while their vertical components are equal by magnitude and opposite. In gen- 
eral, the vectors of the two velocities are perpendicular to each other, implying that the 
mode is hydrodynamically transverse, i.e., the fluid elements oscillate parallel to the lines 
of constant phase. From Eq. ( 125|) it follows that the pure IG mode cannot propagate 
in the direction of the gravity vector (the phase velocity vanishes for = 0), while its 
group velocity vanishes for kz = (no energy transfer by the IG wave for its purely 
horizontal propagation). For an oblique propagation, if the phase velocity has a positive 
2;-component the corresponding group velocity component is negative, and vice versa. In 
the laboratory experiments with stratified liquids, e.g., a water with some vertical salinity 
profile [see Ref. 5 and references cited therein], any vertical or horizontal displacement 
due to some source will cause oblique IG waves, with vertical group velocity component 
dependent on the frequency. It turns out that the more frequent oscillation of the source 
the more vertical transfer of the wave energy. An interesting application of these features 
to the heating of chromosphere by the IG waves has been proposed long ago.^^ Very recent 
observations in fact confirm that model, showing that almost half of the required energy 
for the chromospheric heating comes from the IG waves launched from the photosphere. 

However, its most important feature for the present study follows for a purely hor- 
izontal mode propagation. In such a case the perturbations propagate in the direction 
in which the medium is isotropic, while their z-dependence is self-consistently taken into 
account through Eqs. (I17ti20p . yet the two directions are not mixed, the corresponding 
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characteristic scale-lengths are practically independent and the common exponential term 
then cancels out. In this limit from Eq. it is seen that u = Ub- Hence, the wave 
frequency achieves its maximum value, but more importantly, it becomes independent on 
the wave number. This further implies that, although the wave frequency is constant and 
strictly determined by the temperature and gravity [see Eq. ( 122|) ]. the wave phase speed is 
VjG = ojb/kx- Being decoupled from the lA mode, the horizontally propagating IG waves 
are purely transverse, the actual displacement of the plasma due to the wave is in ver- 
tical direction, and, regardless of their wave-lengths, their frequency is constant. Hence, 
physically the mode is completely different from the fast magneto-acoustic wave which, 
in the presence of a vertical background magnetic field Bq = BqCz, implies the density 
compressions in horizontal direction and the compressive component of the magnetic field 
in vertical ^-direction. 

Exactly because of such a geometry of the horizontally propagating IG mode, and in 
particular because of the vertical displacement of the plasma (i.e., along the magnetic field 
vector), the magnetic field was neglected from the beginning, because a plasma volume 
element moving up and down along the magnetic field lines due to the wave (propagating 
horizontally) will not be affected by the field. Of course, such an idealized picture will 
hold only as long as the frequencies of the two modes are well separated. Otherwise, there 
will be two physically rather different modes (the IG and the fast magneto-acoustic one), 
described by a dispersion equation similar to fl^T]) . but with the magnetic field effects 
included. 

III. APPLICATION TO SLOW EIT WAVES 

To get some feeling about the quantities defined above, adopting Te = = 10^ K, 
7 = 5/3, g = 275 m/s^, one finds Vs = 148 km/s, coc = 1.23 • 10"^ Hz, Ub = 1.07 ■ 10"^ Hz, 
and Hp = 6 ■ 10^ m. Note that the corresponding electron thermal velocity Vre = 
m/s justifies the earlier used inertia-less electron limit. 

Because Ub is so close to Uc, the two modes can only be well separated for large enough 
values of k^. What 'large enough' may mean is seen in Fig. 1, where the ratio of the internal 
gravity and ion acoustic wave phase velocities Vig/via = ^ig/^ia = ^b/{'^1 + k^v'^Y^'^ is 
given in terms of the wave number kx (multiplied by a factor 10'') and for two plasma 
temperatures Tg = Tj = T = 1 and 2 million Kelvin. It is seen that the modes are well 
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separated within a very large span of the wave-lengths that may be of interest for the 
EIT modes. However for wavelengths of the order of 10^ Mm (far left side in Fig. 1) the 
coupling with the acoustic mode should be taken into account. 

For the value of = 4: ■ 10~^ m"^ (i.e., A ^ 150 Mm), from the graph in Fig. 1, 
the phase speed becomes around 25 km/s. This perfectly describes the slow EIT 
perturbations. However, the constant frequency of the mode Ub = 1.07 ■ 10^^ Hz, implies 
perturbations with the wave period T = 2TT/u)b of 98 minutes, that look too long. The 
wave period can become lower for shorter Hp. In fact, the value for Hp used above 
appears to be larger by a factor 2-3 from the values based on observations in coronal 
loops. Hence, taking Hp shorter by a factor 3, the dispersion equation (^T^ becomes 
modified. After solving such a new dispersion equation we find Ub = 0.0026 Hz, and thus 
the wave period T = 2tt /ub = 4:0 minutes. In addition, accounting for the presence of 10 
percent of helium and calculating the mean density, yields the mean ion mass equal to 
1.3 proton mass, and with this the period can further be reduced to 35 minutes. 

Clearly, the buoyancy frequency and the period of the IG wave are sensitive to small 
variations of parameters. The value for the wave period presented here is about the best 
that can be obtained from the present model. An additional improvement would be to 
include the magnetic field and to study the coupled fast magneto-acoustic (FMA) and 
IG modes. The actual displacements of the fluid due to these two modes (and this holds 
for the lA case above too) are in fact orthogonal to each other. So the simple picture of 
the transverse (vertical) incompressible plasma displacement due to the IG mode alone, 
will become much more complicated in the presence of the magnetic field and when the 
two modes are coupled. The coupling to the FMA mode will make the IG mode partly 
compressible too, and in the same time it will affect frequencies of both modes. This 
should particularly hold for a slightly oblique propagation with respect to the vertical 
direction, when the pure IG mode is described by Eq. (125|) . However, the magnetic field 
perturbations cannot be treated exactly within the present model. This is seen from the 
following. Having the magnetic field perturbed implies the use of the Ampere law 

V X 5i = noji = /ioeno (vn - v^i) . (28) 

In view of the earlier introduced variation of the velocity f fT9|) . and for the given equilibrium 
density variation ( [T0|) . this would require that the variation of Bi is the same as the 
variation of the density f fT7|) . ~ exp{—z/2H). However, one is supposed to use the Faraday 
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law too, V X = —dB/dt, where the variation of Bi in the ^-direction should be the 
same as the variation of the electric field potential f l20|) . i.e., ~ exp(2;/2/7). Both are not 
possible, hence there can be no common exponential term to cancel out, and the given 
analytical model does not work. This case can only be studied numerically. 

IV. SURFACE GRAVITY WAVE 

The same procedure for obtaining the slow propagation velocity of the EIT waves 
within the IG wave scheme, can be successfully extended to completely fill in the gap 
between the values obtained above (around 25 km/s in the given example), and the 
minimum possible velocity (i.e., Cg) that would follow from the fast magneto-acoustic 
wave model [which is > 10^ km/s, see Eq. ([T])]. Though this works only for the reduced 
Hp (by a factor 3), and for the mean ion mass equal to 1.3 proton mass, as explained in 
Sec. 3. The maximum wave-lengths are below 200 Mm and wave periods are around 35 
minutes. 

However, even without such modifications, i.e., keeping the original Hp and m^, the gap 
can be filled in also by using the surface gravity (SG) waves in fluids with intermediate 
depths (comparable to the wave-lengths), by identifying the scale length Hp with the 
fluid depth. Performing derivations similar to Sec. 2, the frequency of the gravity wave 
becomes^ 

= [fif/ca; tanh(/ca; iJp)]^/^ 

It should be stressed that a similar model is successfully used also in the description of 
terrestrial Rossby and gravity waves, where the bulk Earth's atmosphere is approximated 
by a layer of the thickness h obtained in the same manner as the quantities H and Hp 
earlier in the text. 

The phase velocity Vsg = ^so/kx of the surface gravity mode is given in Fig. 2 in the 
range of wave-lengths of interest for the EIT waves, for equal to the proton mass, for 
the temperature of the plasma Tg = Tj = 10^ K, and Hp = 6 ■ 10'' m as earlier. Here, the 
wave-length is taken in the interval 0.1 — 160 Mm, yielding the wave periods of the SG 
wave in the interval between 0.8 and 32 minutes, respectively. Hence, the surface gravity 
mode alone is capable of yielding proper propagation speeds for the EIT waves in the 
whole low velocity range. Regarding the coupling with the lA waves, similar arguments 
as in the previous text hold here too. However, because SG wave includes both vertical 
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and horizontal plasma displacement, the coupling with the FMA wave is expected to be 
much more pronounced and it requires numerical simulations. 
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V. CONCLUSIONS 



The present study is focused on the slow range of the propagation speed of the EIT 
perturbations, i.e., those that are well below the magneto- acoustic velocity. Such low 
velocities, that do not fit into the fast magneto-acoustic wave description of the EIT 
perturbations, in the recent past have been a major argument against the wave nature of 
those events. 

In this study, basic plasma and fluid theory is used in order to show that such slow 
propagating perturbations can successfully be described within the theory of gravity 
waves. In fact two possible alternatives are suggested, based on the internal and surface 
gravity wave theory. Hence, these results should at least partly reestablish and support 
the wave-based model of the EIT perturbations. 

In its simplest form, the internal gravity mode described in the text does not produce 
density variation. As already mentioned before, their coupling to the FMA mode in fact 
may include such variation. On the other hand, the horizontally propagating IG mode 
implies vertical displacements (in the z-direction) of the plasma, which means that layers 
of different density (with z) will be mixed together. For example, a denser plasma volume 
raised (within one half-length of the wave propagating in the x-direction) from some level 
Zq to Zi = Zq + 6zq, will be surrounded by a less dense plasma corresponding to the level Zi, 
and vice versa in the other half-length. An observer will see a propagating perturbation 
accompanied by density variations, although the mode itself is not compressible in the 
first place. 

The surface gravity mode is rather different, as it includes both horizontal and vertical 
displacements, and its group and phase velocities are in the same direction. Analysis of 
the fluid motion inside the wave^'^^'^° shows that in the wave crest the fluid motion is 
generally in the direction of the wave, while in the troughs it is opposite. For a horizontal 
propagation the group velocity is non-zero and the mode transfers energy, which is not 
the case with the IG mode. This all should perhaps be a possible recipe for observations, 
to eventually distinguish the SG mode from both the IG and FMA modes. Though as 
a reaction to sudden perturbations (like flares), the coronal plasma can, at any time, 
support different plasma waves that on the other hand may become more or less coupled. 
Lack of ideal conditions, like for example the absence of a clear surface, implies that in a 
realistic situation some hybrid gravity mode may be expected that would include features 
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of both gravity modes discussed in the text. The velocity, and time and space scales of 
the resulting gravity mode should be between the two. A proper description of such waves 
and especially their coupling to FMA mode requires numerical simulations. 
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Figure captions: 



Fig. 1: Ratio of the phase velocities (and frequencies) of the internal gravity mode and 
the ion acoustic mode, in terms of the wave-number. 

Fig. 2: The phase velocity Vsa = '^sc/kx of the surface gravity mode in terms of the 
wave-length. 
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